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ABSTRACT 

The amplitude of density perturbations, for the currently-favoured ACDM cosmol¬ 
ogy, is constrained using the observed properties of galaxy clusters. The catalogue 
used is that of Ikebe et al. (2002). The cluster temperature to mass relation is ob¬ 
tained via iV-body/hydrodynamical simulations including radiative cooling and pre¬ 
heating of cluster gas, which we have previously shown to reproduce well the observed 
temperature-mass relation in the innermost parts of clusters (Thomas et al. 2002). We 
generate and compare mock catalogues via a Monte Carlo method, which allows us to 
constrain the relation between X-ray temperature and luminosity, including its scat¬ 
ter, simultaneously with cosmological parameters. We find a luminosity-temperature 
relation in good agreement with the results of Ikebe et al. (2002), while for the matter 
power spectrum normalization, we find as = 0.78^g'gg at 95 per cent confidence for 
fig = 0.35. Scaling to WMAP’s central value of flo = 0.27 would give a best-fit value 
of as — 0.9. 

Key words: methods: A-body simulations - hydrodynamics - X-rays: galaxies: 
clusters - galaxies: clusters: general 


1 INTRODUCTION 


It has recently become apparent that traditional hydrody- 
namical simulations, where the gas is only allowed to heat 
adiabatically and through shocks, have difficulties in match¬ 
ing observations in the central regions of clusters, with a sig¬ 
nificant underestimation of the temperature corresponding 
to a given cluster mass. This is potentially important for at¬ 
tempts to use the observed temperature function of clusters 
to constrain the matter power spectrum on short scales, a 
topic which has been studied by many authors over the years 
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use hydrodynamical simulations to relate mass to tem¬ 
perature. Such concerns have been g iven further impe¬ 
tus by a recent paper by ISeliakl (' 200211 . who used an ob- 
served relation shi p between cluster temp erature and mass 
(iFinoguenov. Reiprich fc BohrineeilboOlf) . rather than one 


derived from hydrodynamical simulations, to find a normal¬ 
ization for the matter power spectrum signifi cant ly lo wer 
than that of earlier works. In a recent paper llThomas et all 
I 2 OO 2 I 1 we showed that the inclusion of extra gas physics, 
namely radiative cooling of the gas and possible preheat¬ 
ing of the gas before cluster formation, can bring simula¬ 
tions into good agreement with recent Cha ndra observations 
of the cores of clusters llAllen et al . 2J0l, suggesting that 
these may be crucial ingredients in obtaining an accurate 
description of clusters. 

I 11 this paper, we derive a constraint on the matter 
power spectrum normalization as in a way which improves 
on previous work in several ways. On the theoretical side, we 
incorporate the temperature-mass relationship, and its scat¬ 
ter, as obtained from the simulations described above. On 
the observational side, we compare with the data published 
in Ikebe et al. (2002), whose raw catalogue contains around 
one hundred clusters, most with data from both ROSAT 
and ASCA. Finally, on the data analysis side we use a novel 
approach, whereby Monte Carlo simulations are used to gen¬ 
erate mock galaxy cluster catalogues, which through com- 
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parison with the data published in Ikebe et al. (2002) lead to 
a simultaneous constraint on the X-ray temperature to lumi¬ 
nosity relationship, including its scatter, and on the matter 
power spectrum normalization ag. 


2 THE OBSERVED CLUSTER CATALOGUE 

The galaxy cluster catalogue containing the best available X- 
ray data is that compiled by Ikebe et al. (2002) and Reiprich 
& Bohringer (2002). The master catalogue contains 106 clus¬ 
ters, selected by their X-ray ROSAT flux from available clus¬ 
ter catalogues, with 88 among them having been observed 
by ASCA. Imposing a flux cut in the ROSAT [0.1, 2.4] keV 
band of 2.0 x 10 -11 ergs -1 cm -2 , a flux-limited sample of 
63 clusters is then obtained, called HIFLUGCS, which is 
claimed to be statistically complete (Reiprich & Bohringer 
2002). Ikebe et al. (2002) use a slightly different sample in 
their analysis, obtained by excluding the two lowest temper¬ 
ature clusters from HIFLUGCS , ending up with a sample of 
61 clusters with X-ray temperatures ranging from 1.4 keV 
up to 11 keV. Among these, 56 have X-ray temperatures 
derived by Ikebe et al. (2002) from ASCA data by means of 
a two-temperature model that takes into account a possible 
contribution from a cooler component at the cluster core. 

In order to define the observed cluster sample with 
which to compare the artificially-generated cluster cata¬ 
logues, we will impose more restrictive selection criteria on 
HIFLUGCS than Ikebe et al. (2002) did. We will only con¬ 
sider clusters with measured X-ray flux in the [0.1, 2.4] keV 
band above 2.2 x 10 -11 ergs -1 cm -2 , X-ray temperature 
higher than 2 keV, and a redshift between 0.03 and 0.10 
(when performing tests, we found that including clusters 
with z < 0.03 seems to lead to an increase in the best-fit ag 
by a few per cent). These ranges were chosen to maximize 
confidence in completeness of the sample, minimize cosmic 
variance, and because the luminosity-temperature relation 
is expected to deviate from a power law below 2 keV due to 
non-gravitational physics. 

So that we can account for the measurement errors both 
in flux and temperature, which can lead to incompleteness 
effects when imposing either flux or temperature criteria in 
the sample selection procedure, we used Monte Carlo sim¬ 
ulations to generate 40 realizations of the HIFLUGCS cat¬ 
alogue, with the measurement errors in flux and temper¬ 
ature modelled as Gaussian distributed. We then imposed 
our cluster selection criteria, described above, on these cat¬ 
alogues to obtain a set of 40 observed cluster samples, with 
very similar but not identical numbers of clusters, repre¬ 
senting different possible realizations of the chosen observed 
cluster sample. 

We performed extensive tests to determine the mini¬ 
mum number of Monte Carlo realizations of HIFLUGCS 
that should be generated, so as to properly account for the 
effect of the measurement errors on the distribution of the 
cluster properties, within the observed data sample which 
we will use to compare with the artificially-generated clus¬ 
ter catalogues. We found that 40 realizations are enough, 
and increasing their number to 200 or 1000 has a negligi¬ 
ble effect both on the typical distribution of cluster proper¬ 
ties and on the final probability distribution for ag. We also 
generated bootstrap realizations of HIFLUGCS , to deter¬ 


mine whether the flux and temperature measurement errors 
provided by Ikebe et al. (2002) were realistic. They seem 
to be, given that the bootstrap realizations share the same 
mean statistical properties as the Monte Carlo ones, leading 
to negligible differences in the final probability distribution 
of ag values. Finally, there does not seem to be any sys¬ 
tematic shift in the mean statistical properties of both the 
Monte Carlo and Bootstrap catalogues with relation to HI¬ 
FLUGCS , which is reflected by the fact that our result on 
ag does not change even if we just apply our selection crite¬ 
ria to HIFLUGCS, and then compare the resulting cluster 
sample with the artificially-generated cluster catalogues. 


3 THE MOCK CLUSTER CATALOGUES 

The direct simulation of X-ray cluster catalogues from hy- 
drodynamical simulations is beyond present computational 
means due to the excessive number of particles required to 
obtain statistically-robust cluster abundances with temper¬ 
atures above a few keV. Instead, we appeal to the method 
used by Holder, Haiman & Mohr (2001), which is to use 
generalized mass functions of dark matter haloes to gener¬ 
ate catalogues of clusters identified by their redshift and 
mass, and then estimate their X-ray temperatures using 
the mass-temperature relation of clusters in hydrodynam- 
ical simulations. With relation to previous work, the main 
improvement in this paper is the use of a mass-temperature 
relation that is drawn from simulations with more detailed 
models of the intracluster gas physics than have previously 
been implemented, and which closely match the X-ray prop¬ 
erties of observed clusters (Thomas et al. 2002; Muanwong 
et al. 2002, hereafter MTKP02). 


3.1 The mass function 

Currently no standard definition of a dark matter halo ex¬ 
ists, although it is convenient to define a halo as an over- 
dense concentration of matter using the res ults o f the spher¬ 
ical top-hat col lapse model (STH CM; e.g. |PeebletJll99.4 al¬ 
though see also lSheth et a ira . For 0 = 1, the boundary 
of a halo predicted by the STHCM contains a mean internal 
overdensity of 187r 2 « 178 relative to the critical density. 
This result has led many authors to define haloes using an 
overdensity contrast of 200 (which we take as our fiducial 
case). Note that even with the current generation of X-ray 
satellites it is not feasible to measure spatially-resolved prop¬ 
erties of clusters to such large radii. 

A comprehensive study of the mass function of cold dark 
matter haloes was carried out by Jenkins et al. (2001, here¬ 
after JenOl), who compared results from the largest A-body 
simulations available (the Hubble Volume simulations simu¬ 
lated by the Virgo Consortium, which used sufficiently large 
volumes to obtain reliable abundances of haloes on scales 
corresponding to r ich clusters of galaxies) to the mass func¬ 
tion predicted bv [Press fc Schechterl (' 1 974tl . They demon¬ 
strated that the simulated mass function predicts a lower 
abundance of haloes at low masses than the Press-Schechter 
function, but a higher abundance at high masses. Although 
they did not investigate the cause of this discrepancy, they 
pointed out that the Press-Schechter ansatz that all mass is 
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contained in bound objects is untrue in the simulations for 
conventional halo definitions. 

JenOl produced fits to simulated mass functions using 
two different estimators: the splierical-overdensity (SO) and 
friends-of-frien ds (FOF) algorithms. The first case, as imple¬ 
mented bv lLacev fc Colei 1 1994 *), finds and ranks the densest 
dark matter particles and, starting from the densest, grows 
a sphere until the mean internal density equals some multi¬ 
ple of the critical density, p CT , {p ) = Ap cr . Particles within 
this halo are then removed from the list and the procedure 
is repeated until all haloes are fou nd d own t o a given mass 
limit. The FOF algorithm dDavis et alJll985h links particles 
together using a fixed linking length of bn _1 / 3 , where n is 
the mean particle density. FOF does not impose spherical 
symmetry on the shapes of the haloes (which are typically 
triaxial) but can sometimes link together haloes which are 
in close proximity. It is important to use a consistent defini¬ 
tion for cluster masses to define both the mass function and 
the mass-temperature relation: failure to do so can lead to 
errors of 10 per cent in the derived value of ag. 

A further result from the JenOl analysis was that the 
mass functions, when expressed as a function of ln(cr _1 ) 
(where a(AI) is the generalization of ag to any mass-scale), 
are independent of cosmology if haloes are defined using ei¬ 
ther a fixed linking length (e.g. b = 0.2) in the FOF case or 
defining the spherical-overdensity threshold with respect to 
the mean background densit y (e.g . A = 18 012p) in the SO 
case. This was confirmed bv lEvrard et all ( 2002! . hereafter 
Evr02), who also provided fits (as a function of 12) to sim¬ 
ulated mass functions using a SO algorithm with A = 200 
(i.e. overdensity measured with respect to the critical den¬ 
sity). For this paper, we adopt M 200 as the fiducial definition 
of cluster mass, and use the Evr02 fits to estimate the mass 
function at different ( 2 ( 2 ). We have checked that our method 
for measuring cluster masses from the simulations (required 
for the calibration of the mass-temperature relation) pro¬ 
duces almost identical results to the SO method used by 
JenOl and Evr02 (the median difference in halo masses is 
less than 0.5 per cent). 

3.2 The mass—temperature relation 

In this section, we use results drawn from simulations car¬ 
ried out usin g th e hydra 1 V -body/hydrodynam ics code 
dCouchman. Thomas fc Pearcei il995 ; IPearce fc Couchmanl 
Il997t ) on the Cray T3E computer at the Edinburgh Par¬ 
allel Computing Centre as part of the Virgo Consortium 2 
programme of investigations into structure formation in the 
Universe. Details of the method and choice of simulation pa¬ 
rameters were discussed by MTKP02; we summarize details 
pertinent to the results of this paper below. 

We adopt the currently-favoured ACDM cosmological 
model, setting the density parameter I2o = 0.35, cosmolog¬ 
ical constant 12 a = 0.65, baryon density 12b = 0.038, Hub¬ 
ble parameter h = 0.71 and linear power spectrum shape 
parameter F = 0.21. The purpose of this paper is to pro¬ 
vide constraints on ag and so it may seem premature to 
pick one particular value for our simulations. However, the 

1 http://hydra.susx.ac.uk/ 

2 http://virgo.susx.ac.uk/ 


mass-temperature relation of clusters is largely independent 
of ag. The simulations presented in MTKP02 use ag = 0.9; 
we have subsequently repeated one of the simulations with a 
lower normalization, ag = 0.7, and find an identical relation 
within the uncertainties. 

MTKP02 presented 3 simulations which differed in the 
way in which the gas was heated and cooled. In the first 
simulation, a Non-Radiative model, the gas could undergo 
heating by adiabatic compression and shocks but could not 
cool radiatively. Consequently, the resulting clusters are far 
too luminous for their mass and so do not agree with ob¬ 
served X-ray scaling relations (MTKP02). We do not use 
results from this simulation. 

In the Radiative simulation, gas was able to cool radia- 
tive l y us ing the _collisional ionization equilibrium tables of 
[Sutherland fc Dopital ( 1993i f. Cooled material was permit¬ 
ted to form stars, removing low-entropy material with short 
cooling times from the centres of the clusters. Finally, in 
the Preheating simulation (which also includes cooling), the 
specific thermal energy of the gas was raised by 1.5 keV per 
particle at 2 = 4, to crudely model the effects of energy in¬ 
jection by galactic winds. Both models reproduce key X-ray 
cluster scaling relations at 2 = 0, although the former pre¬ 
dicts too much cooled gas (i.e. stars and galaxies) compared 
to observations and the latter too little. 

We estimate the X-ray temperature of each cluster by 
weighting the contribution from each hot gas (T > 10 5 K) 
particle by its bolometric flux 

j, _ TnmjpiAh 0 \(Z, Tj)Tj 

' Y 'Pim.ipiAbo\(Z,Ti) 

Here, mt, pi and Tj are the mass, density and temperature 
of the particles, Z = 0.3 Zq is their metallic ity and Ahni is 
the bolometric cooling function from [Sutherland fc Dopital 
(1993). Adopting a soft-band cooling function (appropriate 
for ROSAT observations) makes no significant difference to 
the estimated temperature. Many clusters show enhanced 
emission from the cluster core that has a lower tempera¬ 
ture than the cluster mean (MTKP02). For this reason, we 
present results for the mass-temperature relation both in¬ 
cluding and excluding the X-ray emission from within the 
‘cooling radius’, defined as the radius within which the mean 
cooling time of the gas is 6 Gyr. The latter results are re¬ 
ferred to as ‘cooling-flow corrected’. 

In Table Q we list parameters for the straight-line rela¬ 
tion of the form 

log(feT/keV) = const + (1/s) log(M 2 oo/h _1 M 0 ) (2) 

that minimizes the dispersion in temperature for all clusters 
with log(Af 2 oo/h _1 M 0 ) > 14. The column labelled “rms” 
gives the root-mean-square dispersion in the log of temper¬ 
ature (for N — 2 degrees of freedom) about the best-fit line. 
We have also measured this dispersion for clusters in a lower 
mass range, 13.7 < log(M 2 oo//i~ 1 M 0 ) < 14, and find very 
similar values. Hence we will assume in our analysis that the 
dispersion is independent of mass. 

The final two columns of Table 1, labelled M 2 oo@3 and 
M 2 00 @6, give the values of the mass, in units of 10 14 /i _1 M 0 , 
for the best-fit relation at temperatures of 3 and 6keV. The 
numbers in the M 2 00 @3 column are mostly very similar to 
each other, except for the top entry for clusters in the Ra¬ 
diative simulation without the cooling-flow correction. The 
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Table 1. Power-law fits to the simulated mass-temperature re¬ 
lations of X-ray clusters: cluster sample; number of clusters in 
sample; slope of relation, s; rms dispersion in temperature about 
best-fit (see text); value of M2oo/10 14 /i —1 Mq at 3keV; value of 
M 2 oo/1O 14 /i- 1 M 0 at 6 keV. 


Sample 



N 

s 

rms 

A/200 @3 

Af200@6 

All data 

Radiative 



36 

1.80 

0.092 

2.9 

10.1 

Preheating , 

&8 

= 0.9 

31 

1.59 

0.056 

2.4 

7.3 

Preheating , 

(78 

= 0.7 

12 

1.75 

0.049 

2.1 

7.1 

Preheating , 

(78 

=any 

43 

1.61 

0.053 

2.4 

7.3 

Cooling-flow corrected 

Radiative 



36 

1.55 

0.079 

2.3 

6.8 

Preheating , 

(78 

= 0.9 

31 

1.51 

0.054 

2.2 

6.2 

Preheating , 

(78 

= 0.7 

12 

1.70 

0.040 

2.4 

7.7 

Preheating , 

(78 

=any 

43 

1.54 

0.049 

2.2 

6.4 


presence of cool gas in the cores of these clusters lowers the 
emission-weighted temperature and hence raises A/ 2 oo@3. 
The slope of the temperature-mass relation for ag = 0.7 
is higher than that for as = 0.9 but the two are in agree¬ 
ment to within the errors; with only 12 clusters covering a 
limited mass-range, the formal 1-sigma error in the slope 
for the as = 0.7 clusters is about ±0.4. The predictions for 
the normalizations of the relations at 6 keV are less certain, 
especially for ag = 0.7, because they require a degree of ex¬ 
trapolation beyond the temperature range of the simulated 
data. For this reason, the difference between the cooling-flow 
corrected normalizations at 6keV for as = 0.7 and as = 0.9 
should not be taken too seriously. We use the combined cata¬ 
logue for our analysis in the next section, but note that very 
similar results are obtained if we use the ag = 0.9 relation 
instead. 

In Tabic |21 we present results from several earlier studies 
of the mass-temperature relation in non-radiative simula¬ 
tions. Note that these results have been obtained by rescal¬ 
ing, when needed, the cluster mass to Af 2 oo (using a NFW 
profile: Navarro, Frenk & White 1995, 1996, 1997) and to 
the cosmology being considered here (as in BN98). Clearly 
there is a wide range of normalizations. This mainly results 
from the different resolutions of the simulations (though in 
the case of EMN96 their method of temperature estima¬ 
tion also plays a part). Also, on average, at fixed temper¬ 
ature the cluster masses in Table [ 2 ] are higher than those 
in Table 1. This is due to the absence of radiative cool¬ 
ing; cluster cores are full of dense, cold gas with short cool¬ 
ing times and this leads to low emission-weighted tempera¬ 
tures. This problem is largely overcome in the Radiative and 
Preheating simulations and can be reduced even further by 
the omission of the cooling-flow component. For compari¬ 
son, Viana & Liddle (1996) and VL99 used a normalization 
for the present-day mass-temperature relation which corre¬ 
sponds to Af 2 oo = 10.1 x 10 14 /i _1 Mg(fcT/6keV) 1 ' 5 , based 
on a simulation of a single higli-mass cluster from White et 
al. (1993b). This agrees well with the results of BN98 and 
ME01, but lies well above the values found by EMN96 and 
in the Radiative and Preheating simulations reported in this 
paper. This change in normalization forces the estimate of 
ag downwards. 


Table 2. Mass-temperature relations_o^2Gra^_clusters from 
previous simulationsj_j5aper_|Evraxd i _Metzlei^^^JiwamJll996|, 
EMN96 ;lBrvan fe Normarll998l BN98: lThoma^^hl200lL T01: 

1; slope of relation, s; value of 
A7 2 oo/10 14 h —1 Mq at 3keV; value of AI 2 oo/10 14 /i~ 4 Mq at 6keV. 


Paper 


s 

M2oo@3 

A/2oo@6 

EMN96 

soft band 

1.50 

2.3 

6.5 

BN98 

bolometric 

1.50 

3.6 

10.2 

T01 

bolometric 

1.50 

2.5 

7.1 

ME01 

bolometric 

1.39 

4.0 

10.6 


3.3 Mock catalogue construction 

We are now in a position to be able to combine the Evr02 
fits to the mass function with the information on the cluster 
mass-temperature relation from the hydrodynamical simu¬ 
lations, to produce mock cluster catalogues with information 
on cluster redshift, mass and X-ray temperature. 

We take the present-day shape of the matter power 
spectrum to be well approximated by that of a cold dark 
matter model, with scale-invariant primordial density per¬ 
turbations and effective shape parameter, T = 0.18. This 
is the favoured value of T from a joint analysis of the 2dF 
(Percival et al. 2001) and SDSS (Szalay et al. 2001; Dodelson 
et al. 2002) data, when accounting for both statistical and 
systematic uncertainties (the allowed interval for T is [0.08, 
0.28] and we confirmed that varying T within this interval 
does not significantly change the final results; changing T to 
either 0.08 or 0.28 leads to a variation of only 0.02 in the 
best-fit ag, with a higher F implying a higher as-). 

We begin by estimating the mean number of clusters 
as a function of mass (M 2 oo) and redshift, using the Evr02 
fits to the mass function, for each value of ag over the inter¬ 
val of interest. Our redshift bins cover the range [0.03,0.10] 
in intervals of 0.001, and our mass bins cover the range 
[0.1, 2.0] x 10 15 A -1 Mq in logarithmically-spaced intervals 
of 0.01. (We have checked that our results are insensitive 
to smaller bin-sizes.) The initial mock cluster catalogues 
are then produced by attributing to each (z, M 200 ) bin, a 
number of clusters drawn from a Poisson distribution whose 
mean is that predicted by the Evr02 fits to the mass func¬ 
tion. We assign a mass and redshift to each individual clus¬ 
ter by randomly drawing the two quantities from a quadratic 
distribution that best reproduces the variation in the cluster 
numbers in the neighbourhood of that bin. In this manner, 
we produce 1000 mock catalogues for each interesting value 
of ag. Through extensive tests we found that such number 
is enough to properly account for the effect of the Poisson 
noise, as increasing the number of mock catalogues per crs 
to e.g. 10000 had a negligible effect on the final probability 
distribution for ag. 

Each cluster in the catalogues is given an X-ray 
temperature, randomly drawn from a Gaussian distribu¬ 
tion in (log 10 A/ 2 oo, log 10 kT), with mean obtained by sub¬ 
stituting the cluster mass [multiplied by H(z)/Ho oc 
-y/T2o(1 ± z) 3 + (1 — flo) to account for the redshift evolu¬ 
tion in the normalization of the cluster X-ray temperature 
to mass relation: see Mathiesen & Evrard 2001] in expres¬ 
sion (2), while the dispersion is assumed to be independent 
of mass. We fix the present-day normalization, slope and 
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dispersion of the mass-temperature relation using the joint 
cluster catalogue obtained from the Preheating simulations, 
where the X-ray temperatures were cooling-flow corrected. 
(Using the parameters deduced from the Radiative simula¬ 
tions does not change the final results significantly.) Our 
method approximately reproduces that used by Ikebe et 
al. (2002) to estimate the observed cluster temperatures. 
We then exclude from the 1000 mock catalogues any cluster 
whose X-ray temperature does not exceed 2 keV. 

To compare our simulated catalogues with the data we 
still need to impose the chosen flux selection criterion, which 
forces us to use a relation between X-ray luminosity (in the 
[0.1, 2.4] keV rest-frame band) and temperature. In order 
to be consistent, we determine this relation from the data 
simultaneously with as (see also Diego et al. 2001). We take 
it to be a power-law of the form 

log 10 (L x //r 2 ergs -1 ) = A + a log 10 (fcT/keV), (3) 

with a dispersion ai ogl0 l x taken to be independent of tem¬ 
perature, and construct a grid of values (with dimensions 
21 x 31 x 16) of the normalization A, slope a and dispersion. 
For each point in this grid, and for every one of the 1000 
catalogues available for each as, we create 50 realizations of 
the luminosity (extensive tests have shown that such num¬ 
ber is enough to lead to a dense coverage of the range of 
possible luminosity distributions, and increasing the num¬ 
ber of realizations to e.g. 200 had a negligible effect on the 
final probability distribution for as) for every cluster by ran¬ 
domly drawing from a Gaussian distribution in (log 10 kT, 
log 10 Lx) with the appropriate mean and dispersion. Every 
cluster then has its X-ray flux in the rest-frame [0.1, 2.4] 
keV band derived, from which the flux in the observed [0.1, 
2.4] keV band is estimated using K-correction formulae. The 
flux limit of 2.2 x 10 -11 ergs -1 cm -2 is then imposed. This 
generates a set of 50000 mock catalogues for each combina¬ 
tion of the four parameters we wish to estimate from the 
data. In all, over twenty five billion mock catalogues were 
generated. 


4 RESULTS 

We are now in possession of an ensemble of catalogues rep¬ 
resenting the observed data set, and a collection of mock 
catalogues for different values of both as and the param¬ 
eters that characterize the X-ray luminosity-temperature 
relation. We chose to perform the comparison between 
the observed and theoretical catalogues via a (three-way) 
2D Kolmogorov-Smirnov test. This test is a generaliza¬ 
tion to two-dimensional distributions of the traditional 
Kolmogorov-Smirnov test, and is due to Fasano & Frances- 
chini (1987), following an earlier idea of Peacock (1983). A 
very good description of the 2D Kolmogorov-Smirnov test 
can be found in Press et al. (1992). In order to calculate the 
probability of each set of 4 free parameters being the correct 
one, we compare each of the observed catalogues with each 
mock catalogue, and then add the probabilities of each pair 
of catalogues being drawn from the same underlying distri¬ 
bution of cluster properties. The probability is taken to be 
zero if the two catalogues being compared do not have the 
same number of clusters, otherwise it is given by the prod¬ 
uct of the probabilities that result from applying the 2D 



Figure 1. Marginalized probability distributions for erg, obtained 
through the 2D Kolmogorov-Smirnov test (full line) and the like¬ 
lihood method (dotted line). 


Kolmogorov-Smirnov test to the 3 available distributions of 
cluster properties — ( z,kT ), ( 2 , Lx) and (kT, Lx). The set 
of free parameters considered most correct will thus be the 
one that most often closely reproduces the observed distri¬ 
bution of the cluster properties (z, kT, Lx). 


4.1 Methodology tests 

As far as we are aware the 2D Kolmogorov-Smirnov test has 
not previously been applied in the same context as here. We 
chose to employ the KS test rather than the widely-used like¬ 
lihood method because it allows a fairly simple incorporation 
of the selection effects entering the observations, and can al¬ 
low for scatter in the cluster relations. Ikebe et al. (2002) 
employed the likelihood function, but did not include scat¬ 
ter in the M-T relation, though Pierpaoli et al. (2003) were 
able to include the scatter in a likelihood analysis. We favour 
the 2D KS test because of its ease of implementation, though 
we do not expect it to lead to significantly different results 
from the likelihood method. 

It is clearly important to compare the 2D KS test to the 
likelihood method. We do this by applying the two methods 
to a simplified situation using 1000 mock observational clus¬ 
ter catalogues (we checked that generating more does not af¬ 
fect the results of the comparison) produced using the Evr02 
fits to the mass function in the same manner as described 
in subsection 3.3, with each cluster being characterized by 
its redshift, z, and mass, A /200 • The assumed fiducial model 
had flo = 0.35, ag = 0.8 and T = 0.18. The sky coverage 
was the same as that of HIFLUGCS and the redshift inter¬ 
val considered was 0.03 < 2 < 0.10. The mock observational 
catalogues were produced assuming that all clusters with 
A /200 above 4.6 x 10 14 fe -1 MQ are detected, and none below. 
In all they have on average 41 clusters (a number similar to 
the HIFLUGCS sub-sample we are working with). 

When applying the 2D Kolmogorov-Smirnov test to 
each of the 1000 mock observational catalogues, 1000 syn¬ 
thetic catalogues were produced for each ag in the interval 
of interest [0.60,1.00], thus overall around 4 x 10' catalogue 
comparisons were made. In this case, the likelihood function 
is the product of the Poisson probabilities of finding exactly 
one cluster in the element dM^oodz at each of the (M^qq, z 1 ) 
combinations present in the mock catalogues, and of find- 
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Figure 2. Marginalized probability distribution for erg, showing 
the bins as calculated and a smoothed version of the distribution. 

ing zero clusters elsewhere in the ( M^oo^z ) plane (see e.g 
Marshall et al. 1983). 

In Fig. 1 we show the probability distributions for erg 
obtained by the two methods. This comparison shows that 
both methods are unbiased, picking up the fiducial as = 0.8 
as the most probable value. Further, the shape of the two 
probability distributions is very similar, though applying 
the 2D Kolmogorov-Smirnov test seems to result in slightly 
more conservative confidence limits. We have made simula¬ 
tions with other initial assumptions and the results do not 
change qualitatively. 

4.2 Application to HIFLUGCS 

The application of the 2D Kolmogorov-Smirnov test to the 
HIFL UGCS data in the manner previously described results 
in the marginalized probability distributions, for each free 
parameter over the three others, presented in Figs. 2 and 
3. The histograms originate from the discretization of our 
parameter space (os, A, a, cri O g 10 L x ) for the Monte Carlo 
simulations. The continuous lines represent the most proba¬ 
ble underlying probability density functions, and result from 
the application of a non-parametric smoothing technique to 
the histogram data. Note that these functions have been 
renormalized for easier comparison with the histograms. In 
summary: 

as ~ 0.78 within [0.72, 1.08] , (4) 

and 

l°gio(Wh -2 ergs -1 ) = A + a log 10 (fcT/keV), (5) 

with 

A ~ 42.1, within [41.2, 42.5] , (6) 

a ~ 2.5 within [1.5, 3.5] , (7) 

ai°g w l x — 0.3 within [0.0, 0.6] . (8) 

where the given ranges are all at the 95 per cent confidence 
level. The most probable combination of the four param¬ 
eters we consider is erg = 0.77, A = 42.2, a = 2.6, and 
°log 10 l x = 0.175. Note that the distribution for os is con¬ 
siderably non-gaussian, with the median value as = 0.81 
being higher than the modal one. The tail extends much 
further to high erg because, as the number of existing clus¬ 
ters increases, it remains possible to reproduce the observed 




a 



Figure 3. Marginalized probability distributions for the normal¬ 
ization A (top) and slope a (middle) of the relation between X-ray 
temperature and luminosity, as well as for its dispersion <Ti ogl0 / x 
(bottom). 


number of clusters by simultaneously choosing lower values 
for A and higher values for ai O g 10 l x ■ In the limit where 
no dispersion in the relation between X-ray luminosity and 
temperature is allowed, the possibility of ag taking high val¬ 
ues disappears, and the marginalized probability distribu¬ 
tion for a g becomes close to gaussian. Re-doing our analysis 
not allowing for any dispersion in the relation between X-ray 
luminosity and temperature, the most probable value for erg 
changes to 0.76, with the 95 per cent confidence interval now 
extending from 0.70 to 0.81, while the most probable values 
for the parameters A and a stay almost the same, changing 
to respectively 42.3 and 2.6. 

These results are similar to those obtained by Ikebe et 
al. (2002). The comparison between the two analysis is made 
difficult by the fact that they only indicate the most prob¬ 
able values for A, a, and <7i O g 10 l x f° r their best-fit flo and 
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as, which are 0.26 and 0.94 respectively if only T > 3 keV 
clusters are considered. Concentrating on this case, and as¬ 
suming Do = 0.26 (plus r = 0.206, as in Ikebe et al. 2002), 
we attempted to recover the values obtained by Ikebe et 
al. (2002) for the other 4 parameters. We took into account 
that their assumed normalization for the cluster mass ( M 200 ) 
to X-ray temperature relation (estimated as if z = 0.05) is 
somewhat higher, such that for a 3 keV cluster they assume 
a cluster mass around 4 per cent higher than we, while at 
6 keV the difference increases to 13 per cent, as well as the 
fact that they do not take into account a possible disper¬ 
sion in the cluster X-ray temperature at fixed mass. A most 
probable value of as = 0.98 was obtained by applying our 
procedure, with good agreement also found for the parame¬ 
ters A, a, and cri ogl0 _l x . Given that some differences remain 
between the two analysis, our results thus seem to be con¬ 
sistent with those of Ikebe et al. (2002). 

In order to determine which type of information in the 
data is driving the results, we determined the most probable 
values for the four parameters under consideration by ap¬ 
plying in isolation the 2D Kolmogorov-Smirnov test to the 
three available distributions of cluster properties — (z, kT), 
(z,Lx) and (kT,Lx)- as is essentially unconstrained by the 
(z, Lx) distribution. All the information comes from the two 
others, with the (kT, Lx) distribution being slightly more 
constraining than the ( z, kT) one. Consistently, the former 
prefers 0.78 as the most probable value for as, while the 
latter settles for 0.79. The information on the parameters 
A, a, and <rio gl0 l x is roughly equally distributed amongst 
the three distributions, though again (kT, Lx) and (z, Lx) 
are always the most and least constraining respectively, and 
when taken in isolation all three distributions lead to very 
similar results. 

In Fig. 4 we compare the cluster properties between a 
realization of the HIFLUGCS sub-sample selected for the 
analysis, the mock sample that most resembles it, generated 
for the most probable set of parameters, and the underly¬ 
ing cluster population. Notice that the incompleteness of the 
flux-limited samples increases considerably as the cluster X- 
ray temperature gets lower, so that below a X-ray temper¬ 
ature of about 5 keV we can conclude that HIFLUGCS is 
vastly incomplete. 

Our analysis is for Do = 0.35, being the value for which 
the large hydrodynamical simulations were run. The spec¬ 
tacular recent results from WMAP (Bennett et al. 2003; 
Spergel et al. 2003) are consistent with this, but their best- 
fit is somewhat lower at Do = 0.27. While we are unable 
to run new large simulations, we can predict the effect on 
ag using the scaling found in earlier analyses. VL99 found 
that for flat cosmologies erg oc D)( 0 ' 47 , and using this scal¬ 
ing we obtain a best-fit as of 0.88 for Do = 0.27. Given the 
small range over which this scaling is needed, the fractional 
uncertainty in erg should be unchanged. 


5 DISCUSSION 

To set the context for the following discussion, we first re¬ 
mind the reader of the constraint from VL99, which for 
Do = 0.35 gave erg = 0.92 within [0.73, 1.12] at 95 per cent 
confidence. By contrast, in Seljak (2002) a value of erg = 0.70 
was obtained, based on the cluster mass to X-ray tempera¬ 
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Figure 4. Comparison of distribution of cluster properties (T vs. 
z, top; Lx vs. T, bottom) between a realization of the HIFL UGCS 
sub-sample considered, the mock sample that most resembles it, 
generated for the most probable set of parameters, and the un¬ 
derlying cluster population. 


ture relation derived in Finoguenov et al. (2001) from cluster 
data. 

The calculation of erg performed in this paper is sub¬ 
stantially different from that carried out in VL99: the semi- 
analytical modelling featured a change in the normalization 
of the assumed cluster mass-temperature relation and in the 
shape of the assumed cluster mass function; a different ob¬ 
servational dataset was used; and the method of estimating 
as from Monte Carlo simulations differed from the previous 
likelihood-type calculation, where no dispersion in the clus¬ 
ter relations was considered. Although, the most probable 
value for erg is quite different in both cases, the 95 per cent 
confidence intervals happen to be very similar. 

In order to find the most important factors behind the 
different results, we ran several Monte Carlo simulations. 
Firstly, we found that the inclusion of dispersion in the 
mass-temperature relation at the level considered in this 
paper does not seem to make much difference. Secondly, 
and more surprisingly, we found that replacing the Evr02 
mass function with the Press-Schechter or the Jenkins et 
al. (2002) mass function also changes the most probable 
value for erg by less than two per cent. This appears in con¬ 
tradiction with claims in the literature, including our own 
(Wu 2001; Pierpaoli et al. 2001; Viana et al. 2002), that the 
choice of mass function can change erg by five to ten percent. 
However that statement is only true if one keeps the mass- 
temperature relationship unchanged, but in fact these differ¬ 
ent mass functions refer to different masses; Evr02 gives the 
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number density of haloes with mass M 200 , Press-Schechter 
uses the virial mass which for the cosmology assumed here 
is about Mi os, and the JenOl mass function corresponds to 
M 63 for the same cosmology. If we use the NFW cluster den¬ 
sity profile to scale these mass functions to the same mass 
definition (e.g. the virial mass), most of the difference in as 
disappears. We note however that this similarity of results 
may be specific to the cosmology adopted here. 

A change that does make a difference is that, as com¬ 
pared to VL99, this paper uses a much lower normalization 
of the cluster mass to X-ray temperature relation. So that we 
could determine the influence of such normalization on our 
results, and be able to compare them more easily with oth¬ 
ers, we calculated the dependence of the most probable as 
on the value of the assumed present-day mass (M 2 oo@ 5 keV) 
of a 5 keV cluster (approximately the median temperature 
of the HIFLUGCS sub-sample we work with). Taking the 
index of the cluster mass to X-ray temperature relation to 
be the standard 1.5, we found 


as = 0.37 + 0.11 x 


/ A/2oo@5keV 


0.83 


(9) 


Note that in our main calculation we assumed M 200 @ 5 keV = 
4.83 x 10 14 /i _ 1 Mq, for an index of 1.54. Given that in VL99 
it was assumed M 2 oo@ 5 keV = 7.67 x 1 O 14 /i _ 1 M 0 , we obtain 
as = 0.97 as the value we would expect from VL99 if the only 
significant difference between the analyses was that change 
in the normalization. Comparing with the VL99 value of 
as = 0.92, this seems to be correct, with the HIFLUGCS 
sub-sample we consider favouring just a slightly higher nor¬ 
malization than the Henry & Arnaud (1991) dataset used in 
VL99. Although it is difficult to untangle all the competing 
effects, we suspect that together the new analysis method 
and the HIFLUGCS dataset allow for a much better esti¬ 
mate of incompleteness which would help explain why they 
favour a higher normalization. 

Turning to comparison with other work, the reason why 
Seljak (2002) obtained a significantly smaller value for as 
with relation to VL99 (to which it is more easily compared), 
is the assumption at fixed cluster temperature of a mass 
that is about 2.4 times lower than that assumed in VL99, 
though this effect is mitigated by Seljak’s assumed local clus¬ 
ter abundance at about 6 keV (from Pierpaoli et al. 2001) 
which was higher than that of VL99. In this paper we too 
have a cluster mass at 6 keV which is much smaller than 
VL99, but the reduction is by a smaller factor of 1.6. 

As we were completing this work, a paper by Pierpaoli 
et al. (2003) appeared in which a similar analysis to ours and 
that in Ikebe et al. (2002) was carried out. The observed clus¬ 
ter sample is also derived from HIFLUGCS , but otherwise 
they use a different approach to obtain constraints on as- 
While both here and in Ikebe et al. (2002) it is attempted 
to constrain as simultaneously with the X-ray temperature 
to luminosity relation, in Pierpaoli et al. (2003) such a rela¬ 
tion is assumed a priori (to be that given by expression 3 in 
Ikebe et al. 2002). We have attempted to reproduce the con¬ 
straint obtained for as by Pierpaoli et al. (2003) when they 
derive the observed cluster sample just from HIFLUGCS. 
Such a constraint can be read from the full line in Figs. 4 
and 5 of Pierpaoli et al. (2003). Concentrating on the case 
of Ho = 0.35, and performing an analysis equivalent to that 
in Pierpaoli et al. (2003), taking care to make the same as¬ 


sumptions and apply the selection criteria in the same man¬ 
ner, but using the 2D-KS method instead, we found a most 
probable value for as and the 90 per cent confidence inter¬ 
val very similar to theirs (a slight overestimation by 0 . 02 ). 
On the other hand, if we just change in our own analy¬ 
sis the cluster mass to X-ray temperature relation so that 
M 2 oo@ 5 keV = 3.82 x 1 O 14 /i - 1 M 0 at present and its index 
to 1.5, as in Pierpaoli et al. (2003), the result for as is again 
very similar (about 0.71) to that obtained in Pierpaoli et 
al. (2003). Clearly, the most significant factor leading to the 
different HIFLUGCS based result obtained here and in Pier¬ 
paoli et al. (2003) regarding as, is the difference in the nor¬ 
malization of the cluster mass to X-ray temperature relation. 

Our results do not indicate a dramatic reduction in as 
derived from the abundance of X-ray clusters. Several other 
recent analyses have favoured low as, for instance from 2dF 
and CMB data (Efstathiou et al. 2002; Lahav et al. 2002), 
from using weak lensing to estimate cluster masses (Viana et 
al. 2002), and from the local X-ray cluster luminosity func¬ 
tion (Allen et al. 2002), but those are at least marginally 
compatible with our present result given the uncertainties. 
Indeed, results from WMAP have forced a modest increase 
in estimates of as via CMB data (Spergel et al. 2003). Our 
estimated value for as is compatible with all published weak 
lensing measurements (e.g Bacon et al. 2002; Hokstra et 
al. 2002; Refregier, Rhodes & Groth 2002; Van Waerbeke 
et al. 2002 ), though only marginally with the very low re¬ 
sults of Brown et al. (2003), Hamana et al. (2002) and Jarvis 
et al. (2003), as well as at the other extreme with that of 
Maoli et al. (2001). 

In the near future, a decrease in the uncertainty in the 
estimation of as from X-ray clusters could come from es¬ 
sentially two sources. On the theoretical side, it would be 
important to reliably estimate the X-ray luminosity of clus¬ 
ters using hydrodynamical TV-body simulations. This would 
enable one to bypass the X-ray temperature as the clus¬ 
ter mass estimator. Though temperature is more reliable, 
it ends up not being as useful as it could be due to the 
fact that all cluster catalogues are flux-limited instead of 
temperature selected, so an estimation of the cluster X-ray 
flux always needs to be made. On the observational side, 
both an improvement in the temperature determination and 
an increase in the range of redshift probed (i.e. a decrease 
in the X-ray flux detection limit) would help bring down 
the uncertainty in the estimation of as- Hopefully, both can 
be achieved with the X-ray satellites Chandra and XMM- 
Newton. In particular, it is expected that the serendipitous 
cluster survey XCS (Romer et al. 2001) to be assembled with 
XMM-Newton data will help greatly in both issues. 
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